Low-temperature structural phase transitions: 
Phonon-like and relaxation order-parameter dynamics 
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A theoretical study on low-temperature structural phase transitions is presented, in which both 
phonon-like and relaxation order-parameter dynamics are contemplated. While the first limiting case 
has been considered previously, the second one is studied here for the first time. Attention is put on 
the low-temperature asymptotics of the temperature dependence of the generalized susceptibility. 
In the relaxation case, it is found ~ (T^ ~T^)~^ for temperature-independent relaxation times in a 
broad region of the temperature-pressure phase diagram. In contrast to the obtained in the phonon- 
like case, this asymptotics is not modified by long-range interactions (dipolar forces, piezoelectric 
effect, etc.). 



I. INTRODUCTION 

Structural phase transitions are usually divided into 
two classes: displacive and order-disorder ones (see, e.g., 
Ref. 1). The labels refer to limiting cases, but still the 
division is convenient. Among other features, it points 
out the order-parameter dynamics: phonon-like in the 
former and relaxational in the latter. By varying a con- 
trol parameter, e.g. pressure or chemical composition, 
the corresponding transition temperature can virtually 
be driven to zero in both cases. Well known examples 
are SrTiOs (displacive)^ and KH2PO4 (order-disorder) 
The theory of low-T displacive phase transition8^i^*S*i has 
been developed since long time ago (see Ref. H for a 
recent review). However, the one concerning to low-T 
order-disorder transitions is still lacking. 

It is worth mentioning the possibility that, depend- 
ing on the temperature regime one considers, the order- 
parameter dynamics of a given transition may evolve 
from one of the above mentioned limiting cases to the 
another one. Let us consider, for instance, a relaxational 
one. With lowering the transition temperature this dy- 
namics may convert into phonon-like because of, e.g., the 
increasing importance of tunnelingiS But quantum effects 
may be important before that such a conversion, if ac- 
tually exists, takes placeiiS Existing theories on low-T 
structural phase transition8^i^*SiL& are therefore clearly 
insufficient to interpret these cases, as far as they are re- 
stricted to phonon-like dynamics, i.e., they neglect any 
damping. The aim of this paper is just to account for 
this damping and, ultimately, to describe low-T struc- 
tural phase transitions with relaxational order-parameter 
dynamics. 

The importance of the order-parameter dynamics in 
low-T phase transitions is due to intimate connection ex- 
isting between statistics and dynamics in quantum sys- 
tems. This has been noticed, for instance, in magnets. 
Here various types of order-parameter dynamics have 
been revealed, and much attention has been paid to the 
role that these dynamics play in the corresponding low-T 
transitions (see Ref. for a recent review) . As we have 



mentioned, different types of order-parameter dynamics 
are also possible in low-T structural phase transitions. 
Here, in contrast, only the phonon-like one has been stud- 
ied until now. In this paper, we shall study both limiting 
cases of phonon-like and relaxation order-parameter dy- 
namics in structural phase transitions, paying attention 
to possible differences in the corresponding anomalies. 
To do so, we shall use a semiphenomenological approach 
in which low-T phase transitions are described within 
a continuous media theory. Within this formalism, we 
shall assume that relaxation is characterized by a single 
phenomenological parameter (a viscosity coefficient) de- 
scribing the simplest case of wavevector and frequency 
independent damping. To determine the microscopic na- 
ture of such a damping is a very difficult task that, being 
the matter of a full microscopic theory, is beyond the 
scope of the present work. So its temperature depen- 
dence, for instance, is assumed to be obtainable from 
the experiments. In the following we shall consider a 
temperature-independent damping for the sake of illus- 
tration, although this dependence would modify final for- 
mulas in a trivial way. 

It is worth mentioning that the lack of studies deal- 
ing with relaxational dynamics in low-T structural phase 
transitions sometimes give rise to some confusion. In 
Ref. il^ for instance, it is reported the low-T thermal 
anomaly of the dielectric constant occurring in prototyp- 
ical order-disorder systems (KH2PO4 and KD2PO4). An 
order-parameter dynamics of relaxation type is well pos- 
sible at these temperatures. Nevertheless, the authors 
analyze the experimental data on the basis of models in 
which this relaxation is completely absent fi^ Their strik- 
ing conclusion that, on the basis of corresponding dis- 
crepancies between theory and experiments, a classical 
order-disorder transition becomes displacive at low tem- 
peratures is therefore unjustified (see below). 

In this work, we focus our attention in the anoma- 
lies of the generalized susceptibility (i.e. the dielectric 
susceptibility when considering ferroelectrics) associated 
with low-T phase transitions because of the current ex- 
perimental interestii^ Nevertheless, it is provided a the- 
oretical framework according to which all other anoma- 
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lies can be obtained. The paper is organized as follows. 
Sec. II has a pedagogical character, and here we con- 
sider a simple model of one-component order-parameter 
in which long-range forces are neglected. It serves us to 
introduce the method to include relaxational dynamics 
when studying low-T structural phase transitions. We 
first reproduce by this method the results of the existing 
(phononic) theory, and then extend this theory to relax- 
ation case. We shall do it in two different ways. One 
of them consists in to calculate the response function 
from the nonlinear equation of motion (Sec. II. A), and the 
other one is by computing the free energy (Sec. II. B). In 
Sec. Ill we improve the model by accounting for acous- 
tic phonons. In Sec. IV we shall focus our attention 
to KDP-type systems. In addition to considering here 
the relaxation case of the order-parameter dynamics, we 
revise the phonon-like case because we found that the 
correct low-temperature asymptotics of the temperature 
dependence of the susceptibility is different from that 
previously reported. 

II. SIMPLE MODEL 

Suppose that the phase transition at zero temperature 
takes place because of change in pressure. Integrating out 
all degrees of freedom but the critical ones, the effective 
potential energy per volume unit of the system can be 
written as 

U-J {lv' + \v'+"^{Vvr)dr (1) 

(it is assumed that the system has unit volume). Here 
r] is the order parameter and a = ap{p — pc) is the only 
pressure dependent coefficient, pc would correspond to 
the phase transition pressure in the classical case, i.e., if 
the ions would have infinite masses and therefore quan- 
tum fluctuations were suppressed. Real phase transition 
pressure is different and will be calculated below. 

Let us here make a comment relative to fluctuations. It 
is frequently said that, while usual (high-T) phase tran- 
sitions are due to thermal fluctuations, that occurring at 
low temperatures are due to quantum ones. The term 
"quantum phase transition" seems then to be synonym 
of low-T phase transition. This makes sense for some 
systems^ but it is quite confusing for others. Among 
these others we have just the structural ones. In fact, as 
we have pointed out, a zero-T structural phase transition 
is perfectly imaginable without any quantum fluctuation: 
a crystal with infinite-mass ions may change its structure 
by applying pressure in accordance with Eq. (^). This 
observation does not mean, as we shall see below, that 
quantum fluctuations and/or quantum statistic play no 
role in low-T structural phase transitions. 



A. Generalized susceptibility 

Let us first calculate the static generalized susceptibil- 
ity following a standard way and considering the classical 
case. The equation of motion for the order parameter can 
be written as 

PcV + jr] + arj + brf — cSJ^-q ~ h, (2) 

where pc is an optic density and 7 a "viscosity coeffi- 
cient." It is assumed that such a viscosity arises because 
of the coupling between the order parameter an other 
degrees of freedom not considered here explicitly. Con- 
sequently 7 may depend, e.g., on temperature. But in 
the following we shall take it as a constant for the sake 
of simplicity. Thermal fiuctuations are described by the 
random field h. In the classical case it is such that 

(/j(r, t)h{v',t')) - 2^frH{v - v')5{t - t'), (3) 

where (3~^ = ksT and (...) denotes statistical average. 

Let us express the order parameter as a Fourier trans- 
form: ri{r,t) = J2kl rfw?7k,aje''' ''e~*"*. Within a first 
approximation in the anharmonicity, the inverse of the 
static susceptibility associated with the macroscopic de- 
gree of freedom (k — 0) can be taken as 

= a + 36^ / duj{ri\^^u,ri--k-^). (4) 
k 

The average in this expression can be computed by using 
the dynamic susceptibility in the harmonic case {b — 0): 



Pc[<jJ + i\i{k)][ijj + i\2{k)] ' 
where w^(/c) = (a + ck"^)/ pc, "1 I Pc and 

Ai,2(fc) = i{7±^[4u;,2(fc)-7']'^'}. (6) 
We then have 

J i^T^r Jo 

Although this analysis is strictly classical, it has been 
shown that the corresponding fully quantum mechanical 
expression follows merely with the rep lacement /cbT — > 
{hu}/2)coih(h(iuj/2) (see, e.g. Refs. [llandil6|). We then 
have 
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X =a+- / -—5- / Im Hfc>^^)coth ?i/3w/2) = a+ — / -— r > ^5 — : ttjti- 

2 r d'^k / 1 V;[l + Ai(fc)/;^] -V^[l + A2(fc)/;y] \ 
^""^ PPcJ i2nf ^ ^[Ai(fc) - ) ' 



where f„ — 27rn/(fi/3), = i^i, and the psi function is 
the logarithmic derivate of the gamma function. [The 
infinite sum in the second form of Eq. appears as a 
result of integration in the complex w-plane.] 



1. Phonon-like dynamics 

If there is no damping (7 — 0) we have Ai = — A2 = 
iujc- The resulting expression in Eq. is 



where 



3hb f k^dk 



Stt'^Pc J u}c{k) 



(9) 



(10) 



and n{u)) = [exp(S/3a;) — is the Bose-Einstein distri- 
bution function. Eq. Q coincides with the correspond- 
ing expression reported in Refs. and0. 

The transition pressure at zero temperature is such 
that a*jj = 0, instead of a = 0. This renormalization of 
the transition pressure is due to zero-point fluctuations. 
The critical frequencies udk) are also renormalizated in 
such a way, but one has to consider higher-order correc- 
tions to obtain it explicitly. Below we shall assume that 
ujcik) represents these renormalized frequencies. Thus it 
is reproduced the vanishing of the optical frequency LodO) 
that takes place within this model as a result of the phase 
transition, i.e., the softening of the optical branch. 

Far from the transition point [ksT <C fnudO)] the in- 
verse of the susceptibility has an exponential dependence 
with temperature oc T^^"^ exp[—hujc{0)/{kBT)]. It is be- 
cause of the "inefficiency" of the thermal activation of the 
optical phonons due to the hardness of this branch. The 
softening of the branch associated with the transition in- 
creases the "efficiency" of this thermal activation. Thus, 
close to the transition [fc^T ^ Hujc{0)] the temperature 
dependence of the inverse of the susceptibility becomes 
ocT^. 



2. Relaxational dynamics 

In the pure relaxation limit (pc —>■ 0) we have Ai(fc) w 
Wc{k) and \2(k) « 7 — vJc{k), where Wc(k) = w^(fc)/7 = 
/C(fc)/7 then describes the inverse of the relaxation times 



in the symmetric phase [see Eq. Q]. In this case, the 
low-T asymptotics of Eq. lO can be obtained by using 
the asymptotic expansion of the psi function (it will be 
valid over the major of the integration interval). Thus 
we get 



(11) 



where 



klh 



khn["f/ujdk)]dk, 



k^dk 



2TrhY J vjl{k)' 



(12) 
(13) 



In the phonon case, the vicinity of the phase transition 
increases the thermal activation of phonons. This leads 
to different behaviors in the temperature dependence of 
the susceptibility depending on the distance from the 
transition point. In the relaxation case, however, the 
distance from the transition point is not so decisive for 
the thermal activation. In consequence, only the 
behavior is obtained in the low-temperature regime of the 
relaxation case for a temperature-independent viscosity 
coefficient. 



Generalized susceptibility revisited: A method 
for further calculations 



Let us now reproduce the above results, but follow- 
ing a different method. The starting point is now the 
calculation of the Landau potential whatever the order- 
parameter dynamics be, i.e. the free energy as a function 
of the critical macroscopic degree of freedom. Once this 
Landau potential is obtained, all the anomalies associ- 
ated with the low-T structural phase transition under 
consideration can be computed. However, as we have 
mentioned, we shall restrict ourselves to calculate here 
the one of the generalized susceptibility. 

Let us now express the order parameter as a Fourier 
transform: 'q{r,t) — X^k ^k(i)e*'^ "". Thus, Eq. iQJ can be 
rewritten showing explicitly its dependence on the above 
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mentioned macroscopic degree of freedom rjo : 



where 



7\ l^kl 



+ br]o ^ f7k?7k'f?-k-k' 

k,k'/0 
b ^ 

+ 1 2^ 'yk^yk'^^k'-'^-k-k'-k", 

k,k',k"7i0 



a 2 & 4 



/C(fc, 770) = a + 35ryQ + cfc^ 



(14) 



(15) 
(16) 



The equations of motion for the Fourier components of 
the order parameter can be written as 



AU degrees of freedom but that associated with rjo are 
integrated out in this potential. However, it describes 
nonequiUbrium states of the system just because, at this 
point, 770 has a nonfixed value. Then one has to reahze 
that what Eq. ifT^ represents is just the Landau poten- 
tial of the system»i 

We are now in a position to reproduce all the anoma- 
lies associated with the phase transition following well- 
known procedures^ Let us consider the anomaly of the 
generalized susceptibility x- In the symmetric phase we 
have = ^" = ip^^/drj^ _^ (hereafter prime de- 
note partial differentiating with respect to 770). Bearing 
in mind that 



(20a) 
(20b) 



?7k + TJyk + Wc(^' '7o)?7k = 0, 



(17) from Eq. (O one finds that 



where w^(A:,?7o) = /C(fc,77o)/pc- Here we have neglected 
the contribution given by the last two sums in Eq. l|T]l. 
Such a neglection is possible by virtue of a weak an- 
harmonicity. But note that we do not neglect this an- 
harmonicity at all: it is the responsible for the depen- 
dence of characteristic frequencies uoc on r/o- The anhar- 
monicity that we are partially neglecting [for instance, 
by omitting the last two terms in Eq. H14|) ] may give 
some contribution to the relaxation which vanishes at 
zero temperature^i Accordingly 7 may depend on tem- 
perature and/or wavevector k. But as before we shall 
take it as a constant. 



Phonon-like dynamics 



47r2pc 



where 



mb 



n[uJcik)] ,^2 



,(fc) 



k^dk 
u>c{k) ' 



(21) 



(22) 



and •n{ui) = [exp(/i/3a;) — 1]"^ is the Bose-Einstein distri- 
bution function. [Summation over wavevectors has been 
replaced by integration: w (27r)~^/ dk.] Eq. 
coincides with the one obtained in Sec.II.A.l. 



Let us present the method assuming first that the 
order-parameter dynamics is phonon-like. If the anhar- 
monicity is weak, the equations of motion for the Fourier 
components of the order parameter can be written as 



fyk + wf(fc,?yo)f?k = 0. 



(18) 



Within this approximation the system reduces to a sys- 
tem of decoupled harmonic oscillators. So its free energy 
can be written asi^ 



$(770) = + i^ln{2sinh[n/3we(fc,?7o)/2]} 

I- k 



1 



hiJcik,r]o) 



k 



-> ln{l-exp[-;i/?c^,(fc,r;o)]}. (19) 
^ k 



2. Relaxational dynamics 

Let us now calculate the low-temperature asymptotics 
of the temperature dependence of the susceptibility when 
the motion of the order parameter includes some relax- 
ation. Similar to the exposed in Sec. II. A, let us consider 
the equations 



jyk + 7'7k + i^lik, ?7o)?7k = 0. 



(23) 



The system is therefore a system of decoupled damped 
harmonic oscillators. Making use of the partition func- 
tion of such damped oscillatorsiS the Landau potential 
can be obtained. It reads 
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$(77o) = UoM {ln^/3[Ai(fc,r7o)A2(fc, m)] ^'^ - lnr[l + \^(k, ^^)|v\ - lnr[l + \^{k,^^)lv\ } , (24) 



where F is the gamma function and Ai.2(fc,?7o) are given 
by Eq. © with ujcik) — > Ludk^riQ). 

The inverse of the susceptibihty in the symmetric 
phase can be calculated from Eq. H24|l by differentiat- 
ing with respect to 770- Noting that Ai,2(fc, '7o) satisfy the 
relations 

Ai(fc) + A2(fc) =7, (25a) 
Xi(k)X2(k)=Lul{k), (25b) 

[Eqs. (|25|l are the relations existing between the coeffi- 
cients and the roots of the algebraic equation — 7A + 



One can see that the expression obtained here by cal- 
culating the Landau potential of the system coincides 
with Eq. ||SJ). Consequently, the low-T asymptotics of 
Eq. I|27|l in the phonon-like and relaxation limiting cases 
match with those given in Sec. II. A. 1 and Sec. II. A. 2 re- 
spectively. 

III. STRUCTURAL NONFERROELECTRIC 
TRANSITIONS: ACCOUNTING FOR ACOUSTIC 
PHONONS 

In any system, there always exists a coupling between 
the order parameter r/ and the strain tensor Uij via the 
term rj'^uu in the potential energy, i.e., the striction effect. 
As a result of this coupling, acoustic phonons give a con- 
tribution to the thermal anomalies associated with low-T 
structural phase transitions which may be significant be- 
cause acoustic phonons are excitations of low energy. So 
let us improve the above studied model by accounting for 
the striction effect. Thus the model can describe quali- 
tatively nonferroelectrics phase transitions. 

Taking into account the striction effect, the potential 
energy can be written as 

(28) 

where a is the striction coefficient, K the bulk modulus 



I 

= (A — Ai)(A — A2) = 0, i.e., the Vieta relations for the 
roots of such a equation (see, e.g., Ref. 2Qj], and bearing 
in mind Eqs. (|20|l . it is easy to see that 



A'i(/c) = A^(fc) =0, (26a) 
'i'(fc) = -A^'(fc) - , ,.,y (26b) 



Then we further get 



(27) 

I 

and V the dilatation of the system {v = uu). We shall 
take into account the anharmonicity within the same ap- 
proximation than in the preceding section. Thus, this 
potential energy can be written in Fourier space as 

U = Uo{vo,e-) + ^^K^i3{k,T]o,e)^i,k^j.-k (29) 

(summation over double indices is implied), where 

rr/ \ a + aeo b , K ^ 

Uoim, e) = ^—^0 + iVo + ye - (30) 

with e being the homogeneous dilatation; 

^ s f a + Sbrig + ae + ck^ iar^kX , . 

^^Ak,V^,e)=^ ^^2fc2j, (31) 

with p being the density of the system and vi the lon- 
gitudinal velocity of sound; and — (r/k,Wk), with Uk 
being the k-Fourier component of a longitudinal acoustic 
displacement. 

The system is now a system of linearly coupled har- 
monic oscillators. Let us consider that relaxation enters 
in the dynamics of the order parameter only. If this re- 
laxation is similar to that considered in the preceding 
section, i.e. if it is characterized by a viscosity coefficient 
7, the partition function of the system can be calculated 
(see Ref. IT^ . From this partition function we can write 
down the Landau potential of the system: 



ih 



2ujl{k) 



^[l + \l{k)/u]-^[l 



:.[Ai(fc)-A2(fc)] 



k'dk. 
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^iVo) = C/o(»7o,e) 



n/3[Ai(fc,?7o)A2(fc,??o)A3(fc,?7o)A4(fc,?7o)]^/^ 



^lnr[l + A,(fc,77o)/;^] 



where the set {Ai} is such that 



Ai + A2 + A3 + A4 = 7, 

A1A2 + A2A3 + A3A4 + A4A1 + A1A3 + A2A4 — U)'^ + uif, 
A1A2A3 + A2A3A4 + A3A4A1 + A1A3A4 = ^ujf , 

A1A2A3A4 = ujjujf - /, 



(32) 



(33a) 
(33b) 
(33c) 
(33d) 



with uj'^{k,r]o) — /Cii(fc, 77o)/pc, ^i{k) — ink, and f{k,rio) = (arjok)^ / (pcp) ■ Let us mention that Eqs. H33fl are the 
Vieta relations for the roots of the equation 



A^ - 7A3 + (c.2 + ..2)^2 _ ^^2^ ^ ^2^2 _ ^ ^ 0^ 



(34) 



The characteristic frequencies of the system are therefore —iXi. 

To our purposes, the value of the macroscopic degree of freedom e can be calculated by minimizing Eq. (|30|l . After 
doing so we find that 



Uoivo) 



« 2 I ^0 4 
2% + 

^11 (fc, Vo) = a + SbiijQ + ck'^, 



(35) 
(36) 



where bo = b - a^/{2K) and 61 = 6 - a^/{6K). 

Let us now calculate the generalized susceptibility of the system in the symmetric phase (770 = 0) following the 
exposed in Sec. ILB. Expressions for the derivates of the A's with respect to 770 a-re then required. From Eqs. I|33(l 
one can easily see that A-(fc) = 0. For X'-{k) one finds that, for instance. 



1 



Ai(fc)-A2(fc) 



2uc{k)iv'^{k) 



f"{k) 



Xlik)+Lof{k) 



(37) 



where it has been taken into account that A3^4(fc) — ±iLUi{k). Similar expressions for Aj 3 ^ik) are also obtainable from 
Eqs. H33I) . Moreover the roots Ai,2(fc) satisfy the relations (|25|1 . Bearing these relationships in mind we find that the 
inverse of the susceptibility, in addition to Eq. pifl .^-^ has now the contribution 



Ax-' 



r(fc) 

47r3/3 7 uj^{k) 



2u;2(fc) 



7/.[l + Al(fc)/7/] 7/>[l + A2(fc)/7/] 

u;f(k) ' iy[\i{k)-X2{k)]\Xl{k)+ujf{k) XUk)+Luf{k) 



2c.2(fc) 



Im 



ip[l + iLUi{k)/iy] 



k^dk. 



In the purely phononic case, 7 = (Ai 

fo2 



vuji {k) V [Ai {k) -iuji{k)] [A2 {k) - iuji{k)] ^ 

— iwc), Eq. H38() reduces to the already known resulli^ 
n[ujc{k)] + l/2 7i^(fc)] + l/2 



A2 



(27r)2p,p J \iu,ik)[io^^{k) - Lofik)] LUi{k)[Lum " ^fik)] 



k^dk. 



(38) 



(39) 



r 



This contribution yields a temperature dependence of the 
inverse of the susceptibility oc T^. This is the leading 
one far from the transition point [/csT ^ hwdO)]. Close 
to the transition [fcsT 3> ftwc(O)], however, the most 
important contribution oc is obtained from Eq. H21|l . 

Substituting the asymptotic expansion of the psi func- 
tion into Eq. (|38|) we get its low-T behavior in the relax- 



ation case (7 7^ 0, ^ 0): 

Ax"' w Aa* + ACT^, (40) 

where Aa* is a constant and AC = — {[a^ / \^hi{K + 
4/i/3)]}C [see Eq. ((TT1) ]. In this case, the importance 
of acoustic phonons is diminished. This is natural be- 
cause, due to the relaxation, the dynamics of the order 
parameter can be activated even far from the transition 
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point. 



where 



IV. PIEZOELECTRIC EFFECT 



Uo 



(43) 



Let us now study low-T phase transitions taking place 
in KDP-type systems. In addition to the striction effect, 
a linear coupling between the order parameter and one 
of the component of the strain tensor is present in these 
systems, i.e., a piezoelectric effect. The role that this 
piezoelectric effect plays in the corresponding transition 
can be revealed considering the potential energy 



U 



A 2 

2"" 



2 

IJ-Uij ) dr 



(V?7)' + ^mv^ 



dUxyT] 



(41) 



where d is the piezoelectric coefficient, ^ the shear mod- 
ulus and \ = K + 2/Lt/3 (for the sake of symplicity, the 
elastic anisotropy has been taken into account partially) . 
Within the same approximation in the anharmonicity 
that in preceding sections, in Fourier space this can be 
written as 



U = Uo[ria,e) + ]-^ ICij (k, ?7o, e)Cj,kC. 



(42) 



with e representing the tensor of homogeneous strains; 



/Coo ~ 0. + ibrjl + aeu + cfc^, 
/Coi = *("%fcx + dky/2), 

= A^ao = iiamh + dkx/2), 
/Co3 = K^lo = iavokz, 



(44a) 
(44b) 
(44c) 
(44d) 

(i,j = 1,2,3), (44e) 



with vt being the transversal velocity of sound; and 
€k = ('^k, Uk), with Uk being the k- Fourier component of 
the displacement vector. Dipolar long-range interactions 
that also take place in these (ferroelectrics) systems have 
been neglected. The possibility of such a neglection is 
explained below. 

Let us assume, as in previous section, that if were no 
coupling between the order parameter and the elastic de- 
grees of freedom {a = d = 0), the relaxation would be 
present in the order-parameter dynamics only. This re- 
laxation being characterized by a viscosity coefficient 7, 
the Landau potential can be written as 



. TT,^S^)^ ^/3[Ai(fc, 7?o)A2(fc, 7?o)A3(fc, 7?o)A4(fc, ?7o)A5(fc, 77o)A6(fc, ??o)]^/^ -A 1 
'^M -Uo + -}_^^ln co^Mk) h ^""^t' + ''"^/"^ / ' 

(45) 

where the set {Xi} corresponds to the roots of the equation 

A^ - 7A'"^ + {ujI + ujf + ujI)X^ - {cjf + ujt)jX^ + [ujl{ujf + ujt) + ujfuj^ + gi + 2gt]\^ 

-ivfcj^-fX + ujlu^fu:^ +g + gicjf + 2gtcj^ = 0, (46) 

I 

with = ICqo/pc, — vik, ujt = Vtk; and As we have seen in previous sections, zero-point fluctu- 

ations yield a renormalization of the coefficient a. But 
^ djvf - Vt) , k^ + dk k ]k k f47a) *° piezoelectric effect (d / 0), this 

^ p^p ^ '° X y) X y, y j renormalized coefficient a* does not vanish at the r = 

^ phase transition: it takes the value a* = — (i^/(4^). In 

gi = [avokxky + d{ki + /c^)/4], (47b) 

consequence, the optical frequencies (or inverse of relax- 
Q,^^ ation times) which enter the above expressions do not 

2gt = (arjok'^ + dtj:ky). (47c) vanish at the transition point. This point, that must 



PcP 



be taken into account in further calculations, has been 



The values of homogeneous strains can be determined overlooked by previous authors. 

by minimizing Eq. 1431) . After doing so we find that ,1 • rm n ;i t ^ ; .• 1 

01 1 — w o i'rom the expression H45(l tor the Landau potential. 



a + d'^/{4:fi) 2 , ^0 4 



the temperature dependence of the susceptibility can ob- 



Uq = -i 77q, (48) tained in a similar way than the exposed in previous sec- 

, 2 tions. After some straightforward but cumbersome cal- 

ICoo — a -I- bi?]Q + ck . (49) culations [which implies differentiating with respect to 
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TjQ and the use of the Vieta relations for the roots of 
Eq. (|46|) ] one finds that, in the relaxation case, the low- 
temperature asymptotics of such a dependence is oc as 
in the relaxation cases studied previously. Mention that 
complicate dispersions arising from dipolar long-range in- 
teractions do not alter this result. It is because the tem- 
perature dependence is obtained irrespective to integra- 
tion over wavevectors [see, e.g., Eq. This justifies 
the omission of terms related with these interactions in 
Eq. (gg). 

The fact that the long-range interactions do not change 
the low-temperature asymptotic of the temperature de- 
pendence of the susceptibility in the relaxation case is 
because of the following. The elementary excitations of 
lowest energy are that associated with the critical degrees 
of freedom due to the corresponding relaxation. It is the 
case whatever the distance from the transition point be, 
the thermal activation being almost independent on this 
distance. 

If there is no relaxation, the elementary excitations 
of lowest energy are, however, the acoustic ones. We 
have already seen in Sec. Ill a similar situation that 
takes place in nonferroelectrics: (longitudinal) acoustic 
phonons yield the main contribution far from the transi- 
tion point. What is different now is that acoustic phonons 
may yield the main contribution even at the transition 
point. It is because, as we have already noticed, there is 
no vanishing optical frequencies due to the piezoelectric 
effect. Estimations of the corresponding gap show that 
in KDP, for example, optical phonons could be frozen 
out up to temperatures of ^ 10 KiSS We shall consider 
in further calculations that this freezing actually takes 
place, what seems to be well possible many other sys- 
tems (the piezoelectric effect in KDP is not especially 
strong) . Then for the low-T asymptotics of the temper- 
ature dependence of the inverse of the susceptibility one 
finds that 

fa' + (TM)'^/^ a'«^{c\a:\)'/\ 
« I (50) 

where a' = a* — a* and 81.2 are constants.-^ Due to the 
combined influence of both striction and piezoelectric ef- 
fects, the ~ T^/^ behavior close to the transition point 



is obtained. These two effect has not been considered 
both together until now, so the above mentioned asymp- 
totics has been overlooked. Mention that it is related to 
the softeing that, due to the piezoelectric effect, there 
exists in an acoustic branch, i.e., vanishing velocity of a 
transversal acoustic wave. Dipolar interactions do not 
modify drastically such a softening, what justifies their 
omission in Eq. H42|) . 

V. CONCLUSIONS 

We have presented a semiphenomenological theory of 
low-T structural phase transitions by including order- 
parameter dynamics of relaxation type. We have re- 
stricted ourselves to the case in which this relaxation does 
not depend on neither the wavevector nor the frequency. 
Within the continuous media theory we use, the relax- 
ation is then characterized by a single phenomenological 
parameter (viscosity coefficient). Both phonon-like and 
relaxation limiting cases can be reproduced by varying 
this parameter. 

We have focused our attention in the low-T asymp- 
totics of the temperature dependence of the generalized 
susceptibility. A dependence ~ (T^ — T^)~^ close to the 
transition point is commonly accepted as an indication of 
a displacive transition, i.e. a transition involving phonon- 
like order-parameter dynamics, although the presence of 
long-range interaction modifies such a asymptotics J: In 
particular, we have found that due to the combined ef- 
fect of striction and piezoelectric effects (KDP case) this 

asymptotics becomes into ~ (T^/^ — Tc^'^)~^. These 
two effects have not been considered both together until 
now, so this asymptotics is overlooked in previous pa- 
pers. For order-disorder transitions (rclaxational dynam- 
ics), we have found that this dependence is ~ (T^ — T^)~^ 
for temperature-independent relaxation times, irrespec- 
tive to the long-range interactions. So, in principle, in 
KDP-type systems low-T structural phase transitions 
with phonon-like and rclaxational order-parameter, i.e. 
displacive and order-disorder respectively, might be dis- 
tinguishable by these asymptotics; although experimen- 
tally it is quite difhcult. Let us stress that the key 
point for understanding the low-T properties is the order- 
parameter dynamics, what could be determined in ferro- 
electrics by measuring e.g. the corresponding dielectric 
losses. 
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